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ABSTRACT 

Operational EW test and evaluation experiments require that the position of the 
aircraft and other moving objects on the range be known precisely as a function of time. 
Terminal Time-Space-Position Information (TSPI) systems involve the range platforms 
interacting at close distances and therefore require precise trajectory information over a 
restricted volume of space. Terminal TPSI systems are used for tactics evaluation and the 
evaluation of simulated weapons firings (e.g., captive-carry hardware-in-the-loop missile 
simulators). Distributed sensor TSPI systems consist of two or more measurement 
sensors located some distance from each other. Each sensor makes a measurement of 
target angle and range. Distributed sensor systems are more complex than single-point 
systems involving multiple hardware installations, complex mathematical computations 
to extract coordinate information, synchronization of multiple measurements and 
calibration of a number of different stations. 

This paper presents a novel distributed sensor TSPI architecture that provides 
precise positioning information of the target relative to a fixed inertial coordinate system. 
The architecture efficiently integrates the information from an inertial navigation system 
(INS), a global positioning system (GPS) and any number of distributed RF sensors 
which may be located onboard a captive-carry aircraft. The significance of this work is 
that by knowing the target's position in a fixed inertial frame of reference (derived from 
the integration process) an evaluation can be made as to the effectiveness of any 
electronic attack or off-board decoys that might have been launched during the field test 
scenario. The induced INS, GPS and sensor noise and the corresponding errors due to the 
integration process are evaluated numerically as a function of the weapon system being 
used. The accuracy in the targeting information is also quantified and compared with the 
true expected values. 
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I. INTRODUCTION 



A. CAPTIVE CARRY ASM EXPERIMENTS 

Currently, the Naval Research Laboratory, Washington, DC conducts captive-carry 
anti-ship missile (ASM) experiments onboard a modified P-3 aircraft. The purpose of the 
experiments is to evaluate the effectiveness of ship defensive systems against various 
types of anti-ship missiles. Currently, the evaluation process involves collecting 
information pertaining to the time and location of certain events which occur during a 
particular test scenario for post-test analysis. One of the goals of the post-test analysis is 
to accurately locate the P-3, the unit(s) involved in the exercise, and any offboard 
countermeasure devices (i.e., CHAFF, RUBBER DUCKS, etc.). By accurately locating 
the players in the scenario, the effectiveness of the countermeasures can be accurately 
determined. 

NRL's P-3 research aircraft is shown in Figure 1.1. It has the capability of carrying 
eight ASM Simulators beneath its main wings. Each ASM Simulator has an active seeker 
installed. Each ASM Simulator has the capability of providing the following information: 

• Target tracking accuracy 

• Seeker guidance type 

• Transmit frequency 

• Pulse width 

• Pulse repetition interval 

• Scan rate 

• Search sector 
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• Beamwidth 

• Transmit Power 

• Antenna Gain 

• Azimuth, Elevation and Range to target 




Figure 1.1: NRL's P-3 RESEARCH AIRCRAFT 



A typical scenario involves the P-3 aircraft equipped with the ASM Simulators and 
the appropriate data recording equipment flying towards a Naval Surface combatant with a 
self protection suite onboard (e.g., AN/SLQ-32(V)4). The P-3 makes several passes at the 
participating unit from several different aspect angles. Usually, the P-3 begins each run at 
a range of twelve nautical miles from the unit. At this range the seekers are turned on to 
begin the search for the target. As the aircraft approaches the ship, eventually, the seeker 
switches from search to track mode. In track mode the seeker locks on to the target in 
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range and azimuth until the P-3 flies past the target or a countermeasure is employed. 
When a countermeasure is employed, the seeker may or may not break lock on the target. 

B. PRINCIPAL CONTRIBUTIONS 

The algorithms presented in this thesis provide a precise time-space-position 
information system (TSPI) based on the inputs from a number of distributed RF sensors, a 
GPS and an INS. The target's position in a fixed inertial frame of reference is derived 
from the integration of the INS/GPS/Radar systems that are located onboard NRL's P-3 
research aircraft. The significance of this thesis is that it provides the necessary tools to 
help NRL perform an evaluation of the shipboard self protection system effectiveness. By 
knowing the target's position as derived from the integration process an evaluation is also 
made as to the effectiveness of any launched decoys during the test scenario. These 
algorithms help the test engineers determine if the seeker broke lock on the target to track 
an offboard decoy (e.g. chaff). 

Also, the relevance of this thesis applies to other applications as well. By knowing 
a target's position relative to a fixed inertial reference frame, the targeting information can 
be sent to any platform and be understood immediately. This type of targeting information 
is ideal for remote firing of weapons that use GPS as an input. 

This thesis develops new algorithms to integrate the INS and GPS along with the 
distributed onboard seekers for determining the target position relative to a fixed inertial 
frame of reference. Mathematical models of the P-3, INS, GPS, and the ASM seekers are 
developed using SIMULINK. These models are used to provide the necessary inputs for 
the distributed sensor TSPI (DSTSPI) algorithms. Evaluation of the algorithms indicates a 
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significant increase in the accuracy of the target position. The results make evident the 
importance of the INS/GPS integration process in knowing the ASM seeker's position in a 
fixed inertial frame of reference. The induced INS, GPS and RF sensor noise and the 
corresponding errors due to the integration process are evaluated numerically as a function 
of the weapon system being used. The accuracy in the targeting information is also 
quantified and compared with the true expected values. 

C. THESIS ORGANIZATION 

Parts of this thesis assume the reader has a background in the areas of radar and 
aeronautics. Thorough discussions are provided in the areas specifically relevant to the 
algorithms presented in this thesis. 

In Chapter II, the background information needed to gain a full appreciation of this 
thesis is discussed. First, the errors associated with the Global Positioning System (GPS) 
and Inertial Navigation System (INS) are discussed along with basic system descriptions. 
Because both GPS and INS reside within fixed inertial frames of reference, a thorough 
discussion of the fixed inertial coordinate systems that are used in this thesis are provided. 
Also, because of the dynamics of the P-3 during the test scenario, a general methodology 
on how to use lever arm corrections is discussed. Finally, this chapter discusses the 
equations used in the modeling of the monopulse and conical scan seekers that are used in 
the simulation. 

In Chapter III, the DSTSPI algorithms are thoroughly developed. Here, a general 
description of the process is given along with an explanation of the SIMULINK block 
diagrams that are used. 
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In Chapter IV, a discussion is provided on the results of the DSTSPI algorithms. 
Several plots are used to aid in the validation and analysis of the results. In addition, a 
further examination is made as to the utility of the algorithm and whether or not it can live 
up to the user’s expectations based on the relevance to the ASM/EW test scenario. 

Chapter V discusses some important considerations involving the hardware and 
software integration process onboard the P-3. Finally, Chapter VI discusses conclusions 
that are made in the process of this analysis. 
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II. BACKROUND INFORMATION 



A. GLOBAL POSITIONING SYSTEM 

The Global Positioning System (GPS) consists of 24 satellites that orbit the earth. 
Of the 24 satellites three are ready service spares. GPS provides locating information to 
an unlimited number of users. Today, GPS is widely used by both the military and the 
civilian sector. More and more civilians are purchasing inexpensive GPS receivers for 
recreational use. The applications of GPS seem limitless. [Ref. l,p. 121] 

GPS is no doubt the most powerful navigational aid available today. There are, 
however, a few drawbacks to the system. First, the positioning information available from 
GPS is updated at a 10 Hz rate for high-end systems and at a 1 Hz rate for most other 
systems. The update rate becomes important when the application involves real time 
processing. In the case where an update rate greater than 10 Hz is required, the 
positioning information from the GPS receiver must be smoothed by employing a Kalman 
filter or some other predictive algorithm. Some higher-end GPS receivers do provide a 
smoothed output that the user is able to sample at higher frequencies. 

There are several factors that determine the GPS accuracy of each position fix. 
The accuracy can vary based on the time of day, the location of the GPS antenna relative 
to the GPS satellites and the type of GPS receiver used. 
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The major sources of error for GPS are summarized as follows: [Ref. l,p 125] 

• Selective availability 

• Clock differences 

• Receiver noise 

• Multipath 

• Atmospheric delays 

• Ephemeris error 

The sources of error are detailed below. 

1. Atmospheric Delay 

Atmospheric delays occur mainly in the ionosphere and the troposphere. In the 
ionosphere, time delays occur because the charged particles interact with the transmitted 
GPS signals. As a result of this interaction, the GPS signal arrives at the receiver with 
some inherent time delay. Since the rate of propagation is nearly the speed of light, any 
delay in the receipt of the signal translates to an increase in the calculated range. 

The equation for the ionospheric delay is given by 

40.3 

A t=— T TEC (2-1) 

cf 

where At is the delay in seconds, f is the frequency in Hz, c is the speed of light in m/s, 
TEC is the total electron content (electrons/m 2 ) along the signal path and 40.3 is an 
empirically derived constant with units (m 3 /s 2 /eIectrons). Since, in the ionosphere, the 
charged particles are typically between 100 to 1000 kilometers above the earth, the path 
length for interaction can be very large. The pseudorange error can be as great as 40 
meters depending upon the time of day, time of year, the solar cycle and the geographic 
location. [Ref. l,p 126] 
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The lower part of the atmosphere can also cause significant delays in the GPS 
signal. This portion of the atmosphere is called the troposphere. These delays are caused 
by the interaction of the GPS signal with water molecules entrained in the atmosphere. 
Here, a 2 meter error can occur in the pseudorange calculation when a satellite is directly 
overhead and a 28 meter error when a satellite is just above the horizon. There are models 
available that can predict the errors from both the ionospheric and tropospheric delays. 
These models improve the performance of the GPS receiver. 

2. Selective Availability 

Selective availability is the status of GPS when the Department of Defense 
degrades the position accuracy to 100 meters. This is accomplished by dithering the 
individual satellite’s clock signal. Dithering is a process by which the clock signals of the 
satellites are encoded such that only military users are able to take advantage of the 
highest accuracy. 

3. Clock Errors 

No matter how precise and accurate a clock is said to be, all clocks are plagued by 
bias and drift errors. Bias is a fixed error associated with the offset of the clock from the 
correct time. Drift is defined as the rate of change of the accuracy of a clock with time. 
Typical clock errors are on the order of a few nanoseconds. A one nanosecond error 
equates to roughly three tenths of a meter error in the pseudorange calculation. 
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4. Ephemeris Error 

Ephemeris errors occur when a satellite’s position is different from the position 
contained within the navigation message that is sent by each satellite. The Ephemeris 
errors are typically on the order of a half meter in the radial direction and three and a half 
meters in the planar direction. 

5. Multipath 

Multipath errors result when the signal does not take a direct path from the satellite 
to the receiver. A signal may bounce off a particular object before reaching the GPS 
antenna. This causes delays in the reception of the signal which can directly translate into 
errors in the pseudorange calculation. Because of the design of current GPS receivers, this 
effect is usually insignificant and is ignored. 

6. Receiver Noise 

All receivers inherently add noise to the signal they receive. A typical receiver has 
a noise value which contributes around 7.5 meters to the error in the pseudorange 
calculation. The overall error associated with receiver noise is wholly dependent upon the 
quality of the receiver being used. 

One method in eliminating most of the errors inherent with GPS is to use 
Differential GPS. Differential GPS (DGPS) involves having another GPS receiver fixed 
to a known location within close proximity of the user’s GPS receiver. This is called a 
differential GPS station. This differential station becomes susceptible to the same errors 
as the user’s GPS receiver. Because the exact location of the differential station is known, 
the differential station can transmit its location along with the corrections for the errors 
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discussed above. This process greatly improves the accuracy of the positioning system. 
An advertised value for position accuracy for a DGPS receiver is around one meter. 
[Ref. 3] 

The DSTSPI algorithms use a DGPS model for the input of the GPS positioning 
information. This model assumes a position accuracy of one meter with an update rate of 
10 Hz. The development of this model is not discussed in this work. Details are found in 
Reference 1, pages 141-152. 

B. INERTIAL NAVIGATION SYSTEM 

The inertial navigation system (INS) detects changes in the total acceleration of a 
body relative to an inertial frame of reference. The primary component of the INS is the 
inertial measuring unit (IMU). The IMU consists of three accelerometers, three rate gyros, 
and two inclinometers. The accelerometers measure the linear, centripetal and 
gravitational acceleration effects on the body. The INS removes the gravitation 
component of the acceleration from the total acceleration. This resultant acceleration is 
integrated once over time to find the velocity of the body and twice to find the position of 
the body. The rate gyros of the IMU measure the angular velocities of the body. The 
inclinometers measure the orientation of the body relative to the inertial frame of 
reference. This orientation is typically denoted as the pitch, roll and yaw. [Ref. 1, p. 136- 
139] 

There are two types of inertial navigation systems. The first type is the gimbaled 
system. Here, the gimbaled IMU is allowed to rotate about four gimbals while a controller 
maintains the IMU platform in a local level orientation towards true North. The relative 
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differences between the gimbals and the platform provide the inertial quantities. The 
problem with the gimbaled systems is that they are usually very heavy and bulky. 
Because of this, the gimbaled system is not very well suited for an aircraft 
implementation. 

The second type of inertial navigation is the strapdown system. Here, the IMU is 
literally strapped down to the aircraft. In this manner the IMU maintains a constant 
orientation to the local coordinate system of the aircraft. Because of this configuration, the 
output of the IMU is given relative to the aircraft’s coordinates. The INS transforms the 
output of the strapdown IMU to accelerations and angular velocities relative to the inertial 
frame. These transformations require substantially more computing power than what is 
needed for the gimbaled system. Because computers are relatively fast and inexpensive, 
the strapdown system is the one generally found on most aircraft. 

The error sources for the INS are mainly due to errors pertaining to the IMU. 
These errors are categorized as follows: bias, cross-axis sensitivity and noise floor. The 
bias error results when the output has a constant offset from the correct value. If there are 
no accelerations on the body, the total acceleration should be zero. Any other value would 
indicate a bias error. A cross-axis sensitivity error occurs when the IMU is not properly 
aligned with the aircraft’s coordinate system. Because it is virtually impossible to achieve 
a perfect alignment of two mechanical systems, a cross-sensitivity error is always present. 
The noise floor for INS prevents measurements below a specified value. It acts as a 
threshold for the output of the INS. The noise floor is caused by the mechanical and 
electrical components of the INS. 
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The DSTSPI algorithms use a strapdown INS model for the input of the pitch, roll 
and yaw of the P-3. This model assumes an accuracy of 0.5 degrees in the determination 
of the pitch, roll and yaw. The development of the INS model is not discussed in this 
work. Details are found in Reference 1, pages 152-157. 

C. COORDINATE SYSTEMS 

Several coordinate systems are used in the development of the DSTSPI algorithms. 
For this reason, a brief discussion is provided on the basic characteristics of these 
coordinate systems. 

1. Earth-Centered-Earth-Fixed 

The Earth-Centered-Earth-Fixed (ECEF) coordinate system is a right-hand 
orthogonal coordinate system that is fixed at the Earth’s origin [0, 0, 0]. Figure 2.1 shows 
the ECEF coordinate system in relation to the Earth. 




Figure 2.1: Earth-Centered-Earth Fixed Coordinates 
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The x-axis is defined as the direction from the origin through the intersection of 0° 
longitude and the equator. The y-axis is defined as the direction from the origin through 
the intersection of 90° east longitude and the equator. The z-axis is defined as the 
direction from the origin through the point defined as true north, or 90° latitude. The 
ECEF coordinate system rotates with the earth. One important advantage of the ECEF 
coordinate system is that all positions that are referenced to ECEF are independent of the 
model used for the irregularities of the Earth’s surface. [Ref. 2, p. 4] 

2. Geodetic Coordinate System 

The geodetic coordinate system is the coordinate system that is most commonly 
used for navigation. Generally, objects which are fixed to the Earth’s surface are located 
in the geodetic coordinate system by designating their latitude, <(>, and longitude, X. The 
latitude and longitude are defined by their relation to the prime meridian and the equator. 
The geodetic latitude, <{>, is defined as the angle from the equatorial plane to the vertical 
direction of a line normal to the reference ellipsoid passing through a location on the 
reference ellipsoid. The geodetic longitude, X, is defined as the angle between the 
reference plane and a plane passing though a location on the reference ellipsoid defined by 
N, where both planes are perpendicular to equatorial plane. The geodetic height, h, is 
defined as the distance from the reference ellipsoid to the location in a direction normal to 
the reference ellipsoid. The reference ellipsoid is a model that defines the shape of the 
Earth. Most modem navigation systems use the World Geodetic System -- 84 model 
(WGS-84). The WGS-84 model gives the model parameters for the shape of the Earth 
along with the local height correction for a particular latitude and longitude. [Ref. 2, p. 5] 
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Figure 2.2 shows a general depiction of the geodetic coordinate system. Here, N 
is the distance from the polar axis to the surface of the reference ellipsoid, h is the height 
correction for the terrain on the Earth’s surface and P defines a position in geodetic 



coordinates. 




Figure 2.2: Geodetic Coordinate System 



3. Tangent Plane Coordinate System 

The tangent plane coordinate system is a right-hand orthogonal system. The origin 
of the tangent plane coordinate system is defined by passing a plane through a point on the 
earth’s surface which is defined by a particular latitude and longitude. The axes, (X,Y,Z), 
of the tangent plane coordinate system are defined as: [Ref. 2, p. 5] 

• X-axis points toward true North 

• Y-axis points toward true East 

• Z-axis points in the downward direction from the point where the tangent plane 
intersects the reference ellipsoid to a point on the polar axis on the earth 
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Figure 2.3 shows a general depiction of the tangent plane coordinate system. 




Figure 2.3: Tangent Plane Coordinate System 



4. Navigation Coordinate System 

The navigation coordinate system is a right-hand orthogonal coordinate system 
whose origin is collocated with the aircraft’s center of gravity. It maintains a local level 
orientation to the reference ellipsoid in the same manner as the tangent plane. If the origin 
of the aircraft is collocated at the origin of the tangent plane, then the navigation 
coordinate system and tangent plane coordinate system are the same (assuming the 
orientations of the axes are the same.) The origin of the navigation coordinate system is 
defined to be the center of the aircraft’s inertial navigation system. The orientation can be 
chosen to be East North Up (ENU), NED, etc.. [Ref. 2, p. 6] 

Figure 2.4 shows the navigation coordinate system with a North East Down (NED) 
orientation. This is the orientation used in the DSTSPI algorithms developed in this 
thesis. It is convenient to choose NED for the orientation of the axes because the 
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coordinate system naturally aligns itself with true North. The difficulty in dealing with 
NED orientation is that altitude is a negative quantity. 




5. Body Fixed Coordinates 

The body fixed coordinate system is the right-hand orthogonal coordinate system 
that is attached to the aircraft and the axes of which are defined by the input axes of the 
inertial navigation system. The body fixed coordinate system rotates with the aircraft 
which implies that the body accelerations of the aircraft are the same as the body fixed 
accelerations if the inertial navigation system is a strapdown system. Figure 2.5 shows a 
depiction of the body fixed coordinate system on the P-3. The orientation of the axes are 
forward, right wing and down. The orientation is defined as a matter of convention. 
Notice that this orientation is the same orientation used for the navigation and tangent 
plane coordinate systems (e.g., forward in the body fixed coordinates is North in the 
tangent plane coordinates.) Keeping the axes orientations the same for these coordinate 
systems results in simpler coordinate transformations. [Ref. 2, p. 8] 
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Figure 2.5: Body Fixed Coordinate System 



D. LEVER ARM CORRECTION 

To accurately locate the seekers for each ASM Simulator relative to a fixed 
coordinate system such as ECEF, the first problem that must be solved is the fact that each 
seeker is not collocated with the GPS antenna. If they were, the seeker’s location would 
simply be the output of the GPS receiver. Since positioning information from the GPS 
receiver is given relative to the GPS antenna in ECEF or geodetic coordinates, natural 
position errors are inherent based upon the displacement between seeker and GPS antenna. 
In addition, because the P-3 is not, in general, locally level with the earth’s surface (the 
reference ellipsoid), additional errors are inherent because of the pitch, roll and yaw of the 
aircraft. To correct for the position errors that result from GPS antenna displacement from 
each seeker and from the pitch, roll and yaw of the aircraft, a lever arm correction m is 
used. [Ref. 3, discussion] 
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The general algorithm for performing a lever arm correction is shown in Figure 
2.6. Three things must be known before a lever arm correction is made. The first is some 
known fixed inertial position that is located somewhere on the P-3. This is usually the 
position of the GPS antenna given in ECEF coordinates, or the position of the INS given 
in navigation coordinates. The fixed inertial coordinates that are chosen are those in 
which the intended local body position is to be referenced (e.g., if the user wants to 



reference all positions to the tangent plane coordinate system, then this is the coordinate 
system that should be chosen). This position is denoted F1 P PT , where the subscript, PT, is 
defined as the local position (local positions refer to those positions located in the body 



fixed coordinate system which are to lever arm corrected to a fixed inertial coordinate 
system.) An example of a local position is the seeker location. The superscript, FI, is 
defined as the fixed inertial coordinate system in which the local position is to be resolved. 




Figure 2.6: Lever Arm Correction Algorithm 

To accurately perform a lever arm correction, the positions of the GPS antenna and 



the individual seekers for each ASM Simulator must be known relative to the center of the 
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aircraft in the body fixed coordinate system. These measurements are extracted from the 
existing drawings of the P-3 or taken manually. The overall accuracy of the lever arm 
correction depends upon the accuracy and precision of the measurements of the relative 
positions of the distributed sensors. Each simulator has two position vectors that must be 
determined in body fixed coordinates. One position vector describes the location of the 
associated distributed seeker and the other describes the after position of the ASM 
simulator. These two position vectors are used later in the analysis to determine the 
attitude of the Simulator. Notice that the forward location of the ASM simulator is 
essentially the location of the seeker. These position vectors are notionally depicted in 
Figure 2.7. 




Once the position vectors for the GPS antenna and the ASM Simulators are known 
in the body fixed coordinate system relative to the origin of the P-3 (which is the location 
of the INS), the next step is to compute the rotation matrix between the body fixed 
coordinate system and the tangent plane coordinate system 1 . This is done by first knowing 
the pitch, roll and yaw of the aircraft. The pitch, roll and yaw of the aircraft are taken 



20 



from the inertial navigation system (INS). The pitch, roll and yaw are defined as the Euler 
angles between the body fixed coordinate system and the tangent plane coordinate system 
when the origins of the two coordinate systems are collocated. We define pitch as 0, yaw 
as ¥, and roll as d>. [Ref. 1, p. 17] 

In order to derive the rotation matrix between the body fixed coordinate system 
and the tangent plane coordinate system the following notation and assumptions are used: 

• {T} represents the inertial tangent plane coordinate system. 

• {B} represents the body fixed coordinate system. 

• The INS is located at the center of gravity of the P-3 aircraft or the origin of 
{B}. Any displacement errors of INS from the true center gravity of the P-3 are 
assumed to be negligible. 

• The {B} coordinate system is defined with X B as the direction forward of 
motion with Y B pointing in the starboard direction (right wing) and Z B pointing 
in the downward direction. 

• The {T} coordinate system is defined with X T in the direction of true north, Y T 
in the direction of east and Zj pointing down towards the center of the earth 

• Roll, d>, is positive when the rotation about X B is clockwise when looking in 
the positive X B direction. 

• Pitch, ©, is positive when the rotation about Y B is clockwise when looking in 
the positive Y B direction. 

• Yaw, 'P, is positive when the rotation about Z B is clockwise when looking in 
the positive Z B direction. 

The Euler angles are depicted in Figure 2.8. Using the right-hand rule defines the 
direction for rotation for roll, pitch and yaw. 



'Note that the tangent plane and navigation coordinates systems differ only by the altitude of the aircraft. 
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Figure 2.8: Directions Euler angles relative to body fixed coordinates 

When finding the rotation matrix, the order in which the rotation is taken is 
important because Euler angles are not vectors. The sum of d>+© is not generally equal to 
0+0. For purposes of this derivation the order is taken to be Z-Y-X, or rotate about Z 
first, then Y and X last. [Ref. 1, p. 18] [Ref. 2, p. 13] 

It is convenient to collocate the origins of the body fixed and tangent plane 
coordinate systems as seen in Figure 2.9. By collocating these coordinate systems, the 
trigonometric relationships between them are easily found. 
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Figure 2.9: Collocated body fixed and tangent plane coordinate systems 

Given a position vector, B P, which is resolved in {B}, the goal is to find the components of 
B P that are resolved in {T}. To do this, first rotate the {B} coordinates about Z B . The 
angle by which {B} must be rotated is defined by T. Therefore, the first rotation matrix is 
a function of Figure 2.10 shows the rotation about the Z B axis. 

Y' b Y b 




X\ 



► x„ 



OUT OF PAGE 

Figure 2.10: {B} rotated about Z B 
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After the rotation through the yaw angle the new axes are defined as (x' B , Y B ,Z B ) . 

The trigonometric relationships between the unprimed and primed frames are used to 
derive the transformation matrix between the two frames. These trigonometric 



relationships are depicted in Figure 2.11. 
- X B sin ¥ 



Y' Y 

1 B y 1 B 




Figure 2.11: Trigonometric relationships when rotated about Z B 



From Figure 2.11, we see that the components of the new primed system are a 
function of the unprimed system and y. X' B is simply the result of the projection of X B 
onto X' B and the Y B sinT term. This results in Equation (2-2). Equation (2-3) is derived in 
a similar fashion. Also, since the rotation occurred about the Z B axis the z components do 
not change during the rotation. This is summarized by the following relationships: 

X B = X B cos T + Y b sin T (2-2) 



Y b =-X b sin T + Y b cosy 



(2-3) 
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(2-4) 



Z'B ~ 



In matrix form equations (2-2), (2-3) and (2-4) simplify to 



X" 




cosT 


sinT 


o” 


X“ 




Y b 


= 


-sinT 


cosT 


0 


Y b 


(2-5) 


Z’b. 




0 


0 


1 


- Z b. 





Now define the rotation matrix for T as 



T(T) = 



cosT sinT 0 
-sinT cosT 0 
0 0 1 



(2-6) 



Next, rotate the primed coordinates about Y’„ through the angle 0, which is the 
pitch angle of the aircraft. This rotation is seen in Figure 2.12. 



X" B X' b 




OUT OF PACE 



Figure 2.12: {B} rotated about Y’ B 



The trigonometric relationships for this rotation are defined in Figure 2.13. 
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- Z’ B sin © 




Z 1 B cos © 

Figure 2.13: Trigonometric relationships when rotated about Y’ B 



From Figure 2.13, the double primed coordinates is defined by the following 
relationships: 



X B = X’ B cos© - Z B sin© 



y;=y b 

Z B = X B sin 0 + Z B cos© 
In matrix form these equations become 

X 

Y” 

Z' B 





cos© 


0 


-sin© 


X" 


= 


0 


1 


0 


y b 




sin© 


0 


cos© 


Zb. 



(2-7) 

( 2 - 8 ) 

(2-9) 



( 2 - 10 ) 



with the rotation matrix for 0 as 
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T(®) = 



cos© 0 -sin© 
0 1 0 
sin© 0 cos© 



( 2 - 11 ) 



Next, rotate the double primed coordinates through the angle O, which is the roll 
of the aircraft. The result is depicted in Figure 2.14. 




Figure 2.14: {B} rotated about X”, 



The trigonometric relationships between the double primed and triple primed coordinates 
are depicted in Figure 2.15. 
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Y" b COS0 



Z" B sin 




Figure 2.15: Trigonometric relationships when rotated about X” B 



From Figure 2.15, the triple primed coordinates are defined in terms of the double primed 
coordinate and the angle O. The result is the following equations: 



x b = x b 



Y b = Y b cos0 + Z B sin0 



Z B = -Y B sin0 + Z B cos0 



In matrix form, these equations become 

'x; 



1 0 0 
0 cos 0 sin 0 

0 -sinO cos 0 



(2-12) 

(2-13) 

(2-14) 



(2-15) 



with the corresponding rotation matrix for O defined as 
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t(<j>) = 



0 



1 0 
0 cosd> sin<J> . 
0 - sin cj> sin d> 



(2-16) 



Combining Equations (2-5), (2-10) and (2-15) and adhering to a [Z-Y-X] rotation 
sequence results in the following equation: 



"x" 




’l 


0 


0 


COS0 


0 


- sin© 


cosT 


sinT 


o" 


X 


Y 


= 


0 


cosO 


sin O 


0 


1 


0 


-sinT 


cosT 


0 


Y 


Z 


{B} 


0 


— sinO 


COSO 


sin© 


0 


cos© 


0 


0 


1 


Z 



For a lever arm correction, the lever arm coordinates are known in the body fixed 
coordinates. Equation (2-17) applies to knowing a position vector in tangent plane 
coordinates and then finding the position in the body fixed coordinates based on the Euler 
angles. To find the lever arm position vector resolved in the tangent plane coordinates, 
find the inverse of Equation (2-17). The matrix form that results is conveniently written as 



X 




1 


0 


0 


cos® 


0 


-sin© 


COST 


sinT 


o’ 




X 


Y 


= INV< 


0 


cosO 


sinO 


0 


1 


0 


- sin T 


cosT 


0 




Y 


Z 


(T) 


0 


-sinO 


cosO 


sin® 


0 


cos® 


1 


0 


1 




Z 



-J(B) 



.(2-18) 



Now, define the rotation matrix, JR as 





1 


0 


0 


cos© 


0 


-sin© 


cosT 


sinT 


of 


^R = INV 


0 


cosO 


sinO 


0 


1 


0 


-sinT 


cosT 


0 




0 


-sind) 


cosO 


sin® 


0 


cos© 


1 


0 


lj 



which is the rotation from the body fixed coordinates to the tangent plane coordinates. 
This is more compactly written as 

l R = INV {t(o)t(©)t(t)} . (2-20) 
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Now, by taking a known position that is resolved in tangent plane coordinates, say 
the GPS antenna position, and adding this position to a position vector that is rotated from 
the body coordinates to the tangent plane coordinate, the result is a position in the body 
fixed coordinate system that is resolved in a fixed inertial frame of reference. To find a 
position resolved in tangent plane coordinates given a position resolved in body fixed 
coordinates with the origins collocated the following equation is applied 

t P 1 =IRx b P 1 (2-21) 

where the subscript 1 means a point somewhere in the body fixed coordinate system. 
Now, by adding this position to the GPS position in the tangent plane coordinate system 
(fixed in the body frame), results in 

T Pi= T P G rs + I RxBp i- (2-22) 

which is the addition of two position vectors resolved in the same coordinate system. 
Equation (2-22) is a simple recipe for performing a lever arm correction. 

E. SEEKER MODELS 

This thesis uses both a conical scan on receive only seeker (COSRO) model and a 
monopulse seeker model to help in the validation process of the DSTSPI algorithms. 
These models are used because they reflect the most common types of seekers found in a 
typical missile inventory throughout the world. The parameters chosen for these models 
are based on typical performance parameters associated with each seeker type. These 
parameters are summarized in Table 2.1 : 
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Estimated Parameter 


COSRO 


MONOPULSE 


Peak Power, P T (kilowatts) 


250 


30 


First side lobe Ratio (dB) 




24 


Frequency, f (GHz) 


9 


17 


Antenna Gain, G (dB) 


33 


23.4 


Half Power Beam Width, 0 3dB (degrees) 


4.5 


8 


Noise Bandwidth, B (MHz) 


10 


10 


Noise Figure, F (dB) 


11 


9.5 


Range resolution (meters) 


100 


20 


Number of Pulses Integrated, n 


100 


100 



Table 2.1: Estimated seeker parameters 



1. Single Pulse Signal to Noise Ratio 



The first step in modeling a seeker is to estimate the signal to noise ratio (S/N) as a 
function of the range between the seeker and the target. First, an assumption must be 
made about the radar cross section of the target, cr. In general, the radar cross section of a 
typical naval surface combatant is of the order of 3000 m 2 . [Ref. 5, p. 129] 

The peak power radiated by any radar transmitter is defined as P T in watts. The 
power density at a given range, R, from an omnidirectional antenna is given by [Ref. 7] 
and [Ref. 4, p. 10] as 



P T e 



-aR 



4ttR 2 



(2-23) 



where a is the attenuation coefficient given in nepers/km. The attenuation coefficient is 
estimated to be 0.055 dB/km at both 9 and 17 GHz. [Ref. 4, p. 278] 
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The power density from a directional antenna with an antenna gain of G x is given by 



P T G T e~ aR 

4rcR 2 



(2-24) 



The power that is intercepted by the target and re-radiated back to the radar antenna is 
given by 



P T G T cre aR 
4nR 2 

where ct is the radar cross section of the target in square meters. [Ref. 4, p. 11] 
The power density received from the re-radiated power is given by 

P T G T oe- 2aR 

(4*R 2 ) 2 



(2-25) 



(2-26) 



since the attenuation and spreading are over twice the distance. 
The signal collected by the antenna is given by 

P T G T A R ge' 2aR 

(4ti) 2 R 4 



(2-27) 



where A R is the effective area of the receiving antenna which is equivalent to 



A = 



471 



(2-28) 



Now, by assuming the transmitted gain of the antenna, G T , is equivalent to the received 
gain of the antenna, G R , the received signal power is approximated by the following: 



P T G T ae 



2 _ -2aR 2 



(4ti) 3 R 4 



', (watts). 



(2-29) 



Now, define the noise power in the receiver as 
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NOISE = kTBF , (watts), (2-30) 

where k is Boltzman's constant, T is the temperature in degrees Kelvin, B is the receiver 
noise bandwidth and F is the receiver noise figure. The temperature is assumed to be 300 
degrees Kelvin. Therefore, kT is 4.14 x 10' 21 watts/Hz at 300 K. The single pulse signal 
to noise ratio is defined as the signal power collected by the antenna divided by the noise 
in the receiver or 

P r G^e- laR X 2 

(4kY R 4 (kTBF) (2 ’ 31) 

Given the single pulse signal to noise ratio as a function of range, the precision in 
the angular measurement for a COSRO and monopulse seeker can be evaluated. 

2. COSRO Seeker 




For a COSRO type seeker, the error caused by thermal noise over the n pulses 
integrated by the servo loop in each coordinate (azimuth and elevation) is [Ref. 4, p. 383] 

®3dB 



<*e = 



k .>/( S/N )m" 



(2-32) 



where k s is the conical-scan error slope, (S/N) m is the single pulse signal to noise ratio, n 
is the number of pulses being integrated and L k is the crossover loss due to beam squint. 
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Figure 2.16: One-way error slope [Ref. 4, p.389] 



The conical scan error slope, k s , and the crossover loss, L k , are shown in Figure 
2.16 as a function of the normalized squint (offset) angle, 0 k /0 3( ,B- It is assumed that the 
beam has a Gaussian shape which results in an optimum error slope of 1.0 at a normalized 
offset angle of 0.5. 

With kj = 1, the error due to thermal noise becomes 



0 



<*6 = 



3dB 






(2-33) 



The scintillation error associated with a conical scan seeker is given by 



a. = 0.225 



3dB 

vf $ k s y 



(2-34) 
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where t^. is the envelope correlation time, f s is the nutation frequency, and p„ is the servo 
bandwidth. Assuming a target with an envelope correlation time t c 0.06 s, a scan 
frequency of f s = 60 Hz and a servo bandwidth P„ = 2 Hz, for an optimized squint angle 



0 k = 0.50 3dB results in a scintillation error of 




(2-35) 



The total error in the precision of the angular measurement for a conical scan 
seeker is found by adding Equations (2-33) and (2-35). The error in the angular 
measurement, 80, or 8c|> is assumed to be Gaussian distributed and is given by 



where a s is a function of the single pulse signal to noise ratio. These equations and the 
parameter estimations found in Table 2.1 are used to develop a model for a COSRO seeker 
type. 

3. Monopulse Seeker 

The next step is to find the precision of the angular measurement for a monopulse 
seeker. The precision of the angular measurement in a thermal noise environment is 



Here, k m is the difference slope and n is the number of pulses being integrated. The value 
for the monopulse error slope is determined from Figure 2.17. It is assumed that the shape 




(2-36) 



a 




(2-37) 
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of the pulse is Gaussian. From Figure 2.17 with a first sidelobe of 24 dB, the value of k m 
is approximately 1.9. 




Figure 2.17: Difference slope versus sidelobe level [Ref. 4,p. 403] 

The error in the angular measurement, 50, or 8<)> is assumed to be Gaussian 
distributed and is given by 



80 = 54 = 





(2-38) 
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where g s is a function of the single pulse signal to noise ratio. These equations and the 
parameter estimations found in Table 2.1 are used to develop a model for a monopulse 
seeker type. 
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III. ALGORITHM DEVELOPMENT 



A. GENERAL DESCIPTION OF ALGORITHM 

This section provides a general description of the DSTSPI algorithms. One of the 
goals of the DSTSPI algorithms is to minimize the errors in the derived target position due 
to the displacement of the GPS antenna from the positions of the distributed ASM seekers 
and the pitch, roll and yaw of the P-3. As mentioned in Chapter II, this is accomplished by 
using lever arm corrections. Another goal of the DSTSPI algorithms is to provide the 
necessary coordinate transformations and corrections to the seeker-to-target information 
resulting in a derived target position in a fixed inertial coordinate system. 

Figure 3.1 depicts the overall functionality for the algorithm. The integration 
algorithms perform the necessary lever arm corrections and compute the pitch and yaw of 
each Simulator. The output from the integration algorithms in conjunction with the target 
data received from the ASM seekers is used to derive the target position in a fixed inertial 
frame of reference. 
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Figure 3.1: General DSTSPI functionality 



B. EXPLANATION OF SIMULINK BLOCK DIAGRAMS 

In the development of the DSTSPI algorithm, SIMULINK® is used as the primary 
tool for the analysis. SIMULINK is a non-linear simulation package developed by The 
MathWorks Inc. for use with MATLAB®. One distinct advantage in using SIMULINK is 
that the results can be viewed in real time or sent to a file for post simulation analysis. 
Another distinct advantage of SIMULINK is that, like functional block diagrams, it is easy 
to see the logical progression of the algorithm. For this reason several SIMULINK block 
diagrams are included in this section. 
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The underlying assumptions for the algorithm are as follows: 



• GPS antenna is located at the [-5m, Om, Om] position in the body fixed 
coordinates and its position is denoted by P ANX . 

• INS is located at the center of gravity of the P-3 which is the origin of the body 
fixed coordinate system and its position is denoted by P CEN . 

• Each Simulator is assumed to be 4 meters in length, which is strictly arbitrary. 

• Simulator #1 is located between the body fixed coordinates [2m, 5m, Om] and 
[-2m, 5m, Om]. The aft position of the Simulator #1 is denoted by P A1 . The 
forward position of the Simulator or the seeker is denoted by P n . 

• Simulator #2 is located between the body fixed coordinates [2m, -5m, Om] and 
[-2m, -5m, Om]. The aft position of the Simulator #2 is denoted by P A2 . The 
forward position of the Simulator or the seeker is denoted by P,^. 

• The target position is denoted by P TAR . 

• The derived target position is denoted by P Tder . 

• DGPS is used with an update frequency of 10 Hz. 



Figure 3.2 shows the relationships of the ASM Simulators and the GPS antenna in the 
body fixed coordinate system. Here, the values for the positions are given in meters. 




F2 



[2 -5 0] 



Figure 3.2: Sensors in the body fixed coordinates (positions given in meters) 
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The notation used in the development of the DSTSPI algorithms is unique and 
needs to be well defined. The following is a list of notational definitions: 

• A superscript in front of the position signifies the coordinate system the 
position is resolved in. For example, B P ANT is the GPS antenna resolved in the 
body fixed coordinates and T P ANT is the GPS antenna position resolved in the 
tangent plane coordinate system. 

• The rotation matrix is denoted as * R where X is the coordinate system from 

which the rotation is occurring and Y is the coordinate system to which the 
rotation is to occur. 

• A primed variable indicates a value with no associated noise added to it. These 
are the actual or true values associated with the flight of the P-3 (e.g., position, 
Euler angles, etc.). 

• An unprimed variable indicates a value with noise added to it. These are the 
derived values (e.g. GPS, INS, etc.). 
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Figure 3.3: DISTRIBUTED SENSOR TSPI Block Diagram 
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Figure 3.3 is the top-level block diagram for the DSTSPI algorithms. Here, the 
“P-3 FLIGHT PROFILER” provides the GPS position at the aircraft’s center of gravity 
( t P CE n)> the INS pitch, roll and yaw(0,C>, v P), the true position of the aircraft ( T Pc EN ) 
and the true pitch roll and yaw of the aircraft (©',0', V F ). Because the GPS position is 
the location of the GPS antenna in tangent plane coordinates, the block “CONDITION 
GPS” converts the center of gravity position ( T P CEN ) to a position located at the GPS 
antenna ( T P ANT ). The “TARGET MOTION” block generates a target in tangent plane 
coordinates ( T Pj AR ) for use in the algorithm and the validation process. The 
“SIMULATOR #1” block is identical to the “SIMULATOR #2” block except for the 
specific values used in the lever arm corrections and the seeker models employed. For this 
reason, the discussion of the DSTSPI algorithms only includes Simulator #1. The 
components of the “SIMULATOR #1” block are discussed in more detail later in this 
chapter. 
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Figure 3.4: P-3 FLIGHT PROFILER Block Diagram 

The “P-3 FLIGHT PROFILER” is depicted in Figure 3.4. The development of this 
part of the model was performed by LT Chuck Gill, et al, as part of a class project in the 
Aeronautical Engineering Department at the Naval Postgraduate School. A thorough 
discussion of this part of the algorithm is found in Reference 6. Although this block is an 
integral part of the analysis, its only function is to provide the outputs required for this 
algorithm which are assumed to be available on the P-3 test flight vehicle. 
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Figure 3.5: P-3 FLIGHT PROFILER CONTROLLER Block Diagram 



The controller used with the P-3 flight profiler required a few modifications to 
make it fully compatible with the DSTSPI algorithm. First, the controller, which is shown 
in Figure 3.5, is modified to allow the aircraft to steer courses in both the starboard and 
port directions. Prior to this modification, the aircraft would only turn in the port direction 
regardless of the ordered course. This correction is made by writing two MATLAB 
programs which corrects the inconsistencies in the flight controller. These programs are 
listed in Appendix A as “HEAD.M” and “OVER2PI.M.” Another problem corrected is 
the instabilities in the controller for steer commands greater than, say, 30 degrees. This 
correction is made by limiting the difference between the ordered steer course and the 
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actual steer course. This gave the controller time to adjust its output to match to required 



response. 




Figure 3.6: TARGET MOTION Block Diagram 

The “TARGET MOTION” block diagram is shown in Figure 3.6. The target 
motion is simulated by assuming an initial displacement from the location of the origin of 
the model’s tangent plane coordinate system (which can be any fixed inertial coordinate 
system). The initial displacement of the target is [5, 0, 0] in nautical miles, with the x 
axis pointing toward true North, the y axis pointing towards true East and the z axis 
pointing downward. The velocity of the target is fixed at 22 knots in a westerly direction 
which translates to -37.131 feet/second in the y direction. By simply multiplying this 
fixed rate by the elapsed time a new target position is calculated for each time step. The 
true position of the target ( T Pj AR ) is then stored in the file “CST.MAT” for post 
simulation analysis. The block is changed to simulate different types of target motion. 
This is done to obtain the results in Chapter IV. 
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Figure 3.7: CONDITION GPS Block Diagram 



The “CONDITION GPS” block diagram shown in Figure 3.7 performs a lever 



to the GPS position at the GPS antenna location ( T P ANT ). As mentioned earlier, the GPS 
antenna is assumed to be at the [-5, 0, 0] location in body fixed coordinates. This is done 
to invoke realism into the algorithm. In a hardware implementation of this algorithm, the 
GPS input would be initially referenced to the antenna location. To perform the lever arm 
correction the INS generated values for the pitch, roll and yaw (©,<!>,¥ ) are used by the 

function "B2NANT.M" to compute the rotation matrix E R . By multiplying the rotation 
matrix and the position B P ANT and, then, adding the result to the GPS position located at 
the center of gravity of the P-3, results in the following: 



arm correction to translate the GPS position at the center of gravity of the aircraft ( T P CEN ) 




(3-1) 
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Figure 3.8: SIMULATOR #1 Block Diagram 



The “SIMULATOR #1” block diagram is shown in Figure 3.8. This block 
diagram performs several functions. First, it performs the integration of INS, GPS and the 
ASM Simulator to derive a target position in the tangent plane coordinate system. Second, 
it performs the same integration again but with true values for the relevant positions and 
Euler angles of the P-3. In this way, the algorithm is validated. Third, this block diagram 
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stores the DSTSPI derived target positioning data and the true derived target positioning 
data for post simulation analysis. 

To derive the target position from DSTSPI, the GPS antenna position ( T P ANT ) and 
the Euler angles from INS are used to find the INS/GPS derived position of 

seeker #1 ( T P FI ) in tangent plane coordinates and the INS/GPS derived yaw (4^ ) and 
pitch (0 SIM ) of the Simulator. This is accomplished by the “P-3 INTEGRATION” 
block as shown in Figure 3.9. Second, this information derived from INS/GPS is then 
used in conjunction with the seeker-to-target information (from the seeker model) to 
derive the target position. The seeker-to-target parameters consist of the range, (R der ), 

azimuth (AZ — v Pg IM ) and elevation (EL -©j IM ) between seeker #1 and the target while 
accounting for the attitude of the Simulator. This information is consistent with the 
typical output of a seeker and is discussed in more detail later in this chapter. 

The true target position ( T Pj AR ), the true P-3 position ( t Pcen) the true Euler 
angles (0’,O , X P ), are used to derive the true seeker position ( T P F , ), the true Simulator 
yaw (% m ) and pitch (0 S | M ) in the “TRUTH MODEL” block. The true position of 
seeker #1, the true yaw and pitch of the Simulator, and the true position of the target are 
used to derive the true seeker-to-target parameters, namely the true range ( R' ), azimuth 
(AZ -4 / g IM ) and elevation (EL-0 SIM ) between seeker #1 and the target. This 
information is used to validate the algorithm in the “TRUE DERIVED TARGET 
POSITION” block. 
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Figure 3.9: P-3 INTEGRATION Block Diagram 



The “P-3 INTERGRATION” block diagram is shown in Figure 3.9. This block 
performs several functions. First, it takes the GPS position referenced to the GPS antenna 
( t P ant ) and lever arm corrects it to the center of gravity of the P-3 ( T P CEN ). Next, it 
takes the GPS position referenced to the center of gravity of the P-3 ( T P CEN ) and lever arm 
corrects this position to the forward and after positions of the ASM Simulator #1. These 
positions are denoted as T P F1 and T P A1 respectively. After these lever arm corrections 
are made, the forward and after positions of Simulator #1 are known in the tangent plane 
coordinate system. Next, the forward and after positions of Simulator #1 are used to 
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compute the yaw ( V F SI>1 ) and pitch (0 SIM )of the Simulator in the tangent plane coordinate 
system. 



Taking the GPS position at the GPS antenna and lever arm correcting it to the GPS 
position at the center of gravity of the P-3, results in the following: 

(3-2) 



T P = T P + T Rx B P 

X CEN i ANT t B ixa x -CEN 



Where, the value jRx P_ CEN > s computed by the function "B2NCEN.M" using the 
instantaneous values for 0, <£>,¥. The negative sign indicates that the lever arm 
correction is taken from the antenna to the center of the aircraft. 

Now, by using the GPS position at the center of gravity of the P-3 resolved in 
tangent plane coordinates and lever arm correcting it to the position at seeker #1, the 
following result is obtained: 



T P = T P + t Rx b P G-3i 

X F1 X CEN^B IV * X F1 \ J J J 

Where the value jRx B P F1 is computed by the function "B2NP2.M" using the 

instantaneous values for 0, <t>, ¥ . 

Likewise, taking the GPS position at the center of gravity of the P-3 and lever arm 
correcting it to the after part of Simulator #1 , results in the following: 

T Pa,= T Pcen+bRx B P a1 (3-4) 

Where the value JRx b P ai is computed by the function "B2NP1.M" using the 
instantaneous values for 0,0,^ . The functions "B2NP2.M" and "B2NP1.M" are found 
in Appendix A. 
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Given the forward and after positions of the Simulator #1 resolved in tangent plane 
coordinates, the yaw and pitch of Simulator #1 is determined. This is accomplished by 
using the “DETERMINE YAW AND PITCH OF SIMULATOR #1” block which is 
depicted in Figure 3.10. 
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Figure 3.10: DETERMINE YAW AND PITCH OF SIMULATOR #1 Block Diagram 



To determine the yaw and pitch of Simulator #1, first, the difference is taken 
between the forward and after positions of the Simulator. This is given by T P F ,- T P A1 . 
This result is a position vector resolved in the tangent plane coordinate system. The length 
of this vector is equivalent to the absolute length of the Simulator. The angles this 
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position vector make with the tangent plane are defined in terms of spherical coordinates 
as shown in Figure 3.11. 

Z 

£ 



o 

Q 




Figure 3.11: Difference between forward and aft Simulator positions resolved in 

spherical coordinates 

From Figure 3.11 we see that the yaw angle is equivalent to the angle 0 SIM and the 
pitch angle is equivalent to© SIM . The pitch angle is defined to be positive from true 
horizontal to a vertical elevation and negative from true horizontal to straight downward. 
Because the tangent plane coordinate system is defined in a NED orientation, the 
Cartesian to spherical coordinate transformation results in the wrong polarity for the pitch. 
This is why the pitch angle is equal to <|> S , M - 90° vice 90° - <(> S i M . The Simulator yaw angle, 
0 S1M is defined in terms of the Eulerian notation, ^ SIM or = 0 SIM . 
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Figure 3.12: TRUTH MODEL Block Diagram 



The block entitled “TRUTH MODEL” is shown in Figure 3.12. This block 
performs the some of the same basic functions as the “P-3 INTEGRATION” block as 
explained earlier in this chapter. First, it takes the true position of the P-3 referenced to 
the center of gravity of the aircraft ( t Pc E n) ar *d lever arm corrects this position to the 
forward and after positions of the ASM Simulator #1. These true positions are denoted as 
t p f , and T P A1 respectively. Next, it takes forward and aft positions of Simulator #1 and 
computes the yaw ( V P SIM ) and pitch (® SIM ) of the Simulator in the tangent plane 
coordinates by way of the “GENERATE TRUE AZIMUTH ELEVATION RANGE 
BETWEEN SEEKER AND TARGET” block as shown in Figure 3.13. Also, the 
difference is taken between the true seeker position and the true target position which 
results in the value T P TAR - T P F - This value is used in Figure 3.13 to compute the true 
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seeker-to-target parameters namely the range ( R ), the azimuth ( AZ - 4^,^ ) and the 
elevation (EL' - © SIM ). 
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Figure 3.13: GENERATE TRUE AZIMUTH ELEVATION AND RANGE 
BETWEEN SEEKER #1 AND TARGET Block Diagram 



Now, by using the true position at the center of gravity of the P-3 resolved in tangent 
plane coordinates and lever arm correcting it to the position at seeker #1, the following 
result is obtained: 
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(3-5) 




Where the value jR x B P’, is computed by the function "B2NP2.M" using the 

instantaneous values for © ,<f> .'F . Likewise, taking the true position at the center of 
gravity of the P-3 and lever arm correcting it to the after part of Simulator #1, results in the 
following : 



Where the value Jr'x b P m is computed by the function "B2NP1.M" using the 

instantaneous values for 0 ,0 .VP . Next, the true values of the forward and after 
positions of the Simulator #1 and the position vector between seeker #1 and the target are 
used in the SIMULINK block diagram pictured in Figure 3.13 to find the true yaw and 
pitch of the Simulator. To determine the true yaw and pitch of Simulator #1, the 
difference between the forward and after positions of the Simulator, T Pp, - T P^, , is used. 
This is a position vector resolved in the tangent plane coordinates. Again, the angles this 
position vector make with the tangent plane are defined in terms of spherical coordinates. 
(Refer to Figure 3.1 1.) Using the same argument that was used in the derivation of the 
INS/GPS derived yaw and pitch, the true yaw angle of Simulator #1 is equivalent to the 
angle VPs, M and the true pitch angle of Simulator #1 is equivalent to ©s, M . These values 
are obtained form the following equations: 




(3-6) 




(3-7) 



®SIM — ^SIM 90 



(3-8) 
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Next, calculate the seeker-to-target parameters by using the value for the difference 
between the true seeker position and the true target position, T P^ AR - T P n , and the true yaw 
and pitch of Simulator #1. The value for t Ptar -T Pfi is a position vector resolved in 
tangent plane coordinates. The angles this vector makes with the tangent plane are defined 
by a Cartesian to spherical coordinate transformation. This is seen in Figure 3.14. 



o 

Q 




- 90 ° 



Figure 3.14: Difference between true seeker and target positions resolved in 

spherical coordinates 

Here, the angle AZ is the true azimuth angle between seeker #1 and the target without 
consideration for the yaw angle of the Simulator. The angle EL' is the elevation angle 
between seeker #1 and the target without consideration for the pitch angle of the Simulator 
and is defined with respect to horizontal. The values for AZ’ and EL are given by the 
following equations: 
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( 3 - 9 ) 



~ ® TAR 

EL = <l>x AR - 90° (3-10) 

To find the true azimuth angle between the seeker and the target, the azimuth must 
be compensated by an amount equal to the yaw angle of the Simulator. The goal is to get 
a value for the azimuth angle that is an expected output of the Simulator. A typical radar 
would not account for its orientation with its surroundings when giving its output. The 
result after compensation is AZ - . The true yaw angle of the Simulator is 

subtracted from the true azimuth angle to maintain consistency with true North. If, for 
example, the Simulator is pointing due East or 90° and the azimuth angle between the 
seeker and the target is also taken to be due East or 90°, then we would expect the 
Simulator to return an angle of zero for the azimuth. 

Now, using the same reasoning, calculate a compensated elevation angle between 
the seeker and the target. This is given as EL’ - ©s, M . 

The true range between the seeker and the target is a easily obtained from the 
Cartesian to spherical coordinate transformation on the position vector T P TAR - T P F1 since 

the range does not care about what the Simulator is doing. This is given as R . 

To validate the algorithm, the “TRUE DERIVED TARGET POSITION” block 
diagram as pictured in Figure 3.15 is used. The goal of this block diagram is to find the 
target position in tangent plane coordinates using true values for inputs. This is 
accomplished be first taking the true seeker #1 to target parameters and correcting these 
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values for the true yaw and pitch of the Simulator. This allows for the extraction of the 
true position vector between the seeker and the target via a coordinate transformation. 




TRUE POSITION 
OF SEEKER #1 



Figure 3.15: TRUE DERIVED TARGET POSITION Block Diagram 



Given the true seeker to target parameters, R' , AZ’ - ¥s, M and EL' - © SIM , first 

these values must be corrected for the true yaw and pitch of the Simulator. This is done by 
adding the respective values in the following manner: 

AZ = AZ - 'Pj IM + ¥ sim (3-11) 

EL = EL - ©sim + ©sim (3-12) 

If this correction were not made, the spherical to Cartesian coordinate transformation 
would result in severe displacement errors from the yaw and pitch of the Simulator. The 
elevation angle, as mentioned earlier, is defined with respect to horizontal. To put the 
elevation term in spherical coordinates simply add 90 degrees, or <j> TAR = 90° + EL . 
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Next, perform a spherical to Cartesian coordinate transformation. The result in 
tangent plane coordinates is the position vector T P TAR - T P F1 . By adding the true position 
of seeker #1 to this quantity, the result is the true target position in tangent plane 
coordinates. This is seen in the following equation: 

T Ptar= T P;ar- T Pf,+ T P;, (3-13) 

The next step in the algorithm development is to talk about the seeker model for 
the conical scan radar. The block diagram entitled “SEEKER MODEL FOR CONICAL 
SCAN RADAR” is shown in Figure 3.16. The three parameters that are relevant to this 
algorithm are the range, elevation and azimuth angle between the modeled seeker and the 
target while accounting for the ASM Simulator's attitude. The output of a real seeker, as 
mentioned earlier, does not account for its surroundings. It only knows the target in a 
relative sense. 
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Figure 3.16: SEEKER MODEL FOR CONICAL SCAN RADAR Block Diagram 



Because the information received from the seeker is not referenced to a particular 
coordinate system, true values for the position of the seeker are used to derive the seeker 
to target parameters. Also, because the attitude of the seeker is not accounted for in its 
output, true values for the yaw and pitch of the seeker are used. 

Referring to Figure 3.16, first the difference is taken between the true seeker 
position and the true target position which results in the value T P TAR - T P n . This value is 
used to help compute the modeled seeker to target parameters. The value for T Pt A r -T Pfi 
is a position vector resolved in tangent plane coordinates. By performing a Cartesian to 
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spherical coordinate transformation on this position vector, the values R’, <J> XAR , and 
0 TA r are obtained. Note, that these values are based on the true positions of the seeker 
and the target as denoted by the prime. 

Next, the value for R is quantized in 100 meter segments to simulate the range 
resolution of the modeled seeker. The output of the quantizer is R GATE , which is the range 

between the seeker and the target that has an error component based on the range 
resolution of the seeker. This quantity is defined as 

^DER = ^GATE (3*14) 

where “DER” is defined as a derived value. 

The MATLAB function "CONTHETA.M" adjusts the precision in the angular 
resolution of the seeker based on the signal to noise ratio and the scintillation as discussed 
in Chapter II. A similar MATLAB function, “THETAD.M,” performs the same function 
as “CONTHEATA.M” but for the monopulse seeker model found in the block diagram 
entitled “SIMULATOR #2” The signal to noise ratio changes as a function of the range 
between the seeker and the target. The angular resolution is normally distributed resulting 
in the value 50 . The value for 50 contributes errors in both directions of the solid angle 

subtended by the radar beam and, therefore, must be added to both 0!p AR and 4 >tar • 

The angle AZ is the azimuth angle between seeker #1 and the target without 
considering the yaw angle of the Simulator. Likewise, the angle EL is the elevation angle 
between seeker #1 and the target without considering the pitch angle of the Simulator and 
is defined with respect to horizontal. The values for AZ and EL are given by the 
following equations: 
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(3-15) 



AZ = e TAR + 50 
EL = tf T AR +59-90° (3-16) 

The azimuth angle between the seeker and the target after compensating for the 
true yaw of the Simulator is AZ - . Again, the true yaw angle of the Simulator is 

subtracted from the azimuth angle to maintain consistency with true North. The elevation 
angle between the seeker and the target after compensating for the pitch angle of the 
seeker is given by EL - ©s, M . 

The final step in the algorithm is derive the target position from the INS/GPS 
derived parameters and the seeker to target parameters. Shown in Figure 3.17 is the block 
diagram entitled “IN S/GPS/RAD AR DERIVED TARGET POSITION”. This block 
performs the same basic functions as Figure 3.15 with a few minor modifications. 

Given the derived seeker-to -target parameters, R der , AZ - ^sim an d ®sim ’ 
these values must be corrected for the INS/GPS derived yaw and pitch of the Simulator. 
This is done by adding the respective values in the following manner: 

^^DER = ~ ^SIM ^SIM = ®DER (3*17) 

EL der =EL-0 sim +0 s1m (3-18) 

Notice that yaw and pitch terms of the Simulator do not fully cancel. The residual values 
are due to the errors involved with INS, GPS and the seeker. These errors translate 
directly into displacement errors in fixing the target’s position in tangent plane 
coordinates. These errors are discussed in more detail in the following chapter entitled 
“VALIDATION AND ANALYSIS.” 
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Figure 3.17: INS/GPS/RADAR #1 DERIVED TARGET POSITION Block Diagram 



The elevation term, EL der , is not in the proper form. To put the elevation term in 
the form of spherical coordinates simply add 90 degrees, or <|> DER = 90° + EL der . This 
results in all values in terms of spherical coordinates. 

Next, perform a spherical to Cartesian coordinate transformation. The result in 
tangent plane coordinates is the position vector T P Tder - T PFi • To derive the target position 
in tangent plane coordinates, add the INS/GPS derived position of seeker #1 to this 
quantity. This is seen in the following equation: 

T P = T P - T P’ + T P (3-19") 

r TDER r Tder r Fl T 17 FI V-' 1 ') 
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Note, the true and derived position vectors of seeker #1 do not cancel. This is due 
to the fact that the derived seeker-to-target parameters are based on the true seeker 
position. 

Next, assume the target has an upper limit on its speed. This means the x, y and z 
components of the DSTSPI derived target position can be rate limited to a specified value. 
In other words, the target position is not allowed to change by more than a specified 
amount. This rate of change is limited to a value of 50 KTS. The z component of the 
target position could have been limited to values close to the surface, but, this would have 
possibly resulted in the miss detection of any airborne decoys such as a chaff cloud. 
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IV. VALIDATION AND ANALYSIS 



All the figures in this chapter were generated using MATLAB. Also, all the 
figures were generated under the following conditions: 

• The P-3’s initial course is 000° true 

• The P-3’s initial altitude is 2000 feet 

• The P-3 ’s initial speed is 500 ft/sec or 296.242 kts 

• The P-3’s initial position in the tangent plane coordinate system is [0, 0,-2000] 
in feet 

• The Target’s initial course in Chapter IV, Section’s A through C is 270° true 

• The Target’s speed is 22 kts throughout the simulation 

• The Target’s initial location in the tangent plane coordinate system is [5, 0, 0] 
in nautical miles 

• All data collection periods were 60 seconds 

A. INS AND GPS ERRORS 

This section looks at the modeled errors from both INS and GPS. Plots of these 
errors were generated using MATLAB. As discussed in Chapter II, there are several 
factors which contribute to the errors seen from INS and GPS. One of the limiting factors 
in the algorithm is the fact that these errors are always present. The biggest contributors to 
the displacement error in the derived target position are due to the errors from INS and the 
seekers. Although GPS errors are important, an error in the seeker's position only results 
in the same order magnitude error in the derived target’s position. Simply stated, because 
the arc length is equal to the angle in radians multiplied by the radius (or the distance 
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between target and the seeker), any angular errors prove to be detrimental to the algorithm. 
These angular errors result from both the seeker’s ability to measure the azimuth and 
elevation between the Simulator and the target and the INS’s ability to measure the Euler 
angles. As a hypothetical case, assume that a particular seeker's position is known to 
within one meter and the pitch, roll and yaw of the Simulator are known to within 0.5°, or 
0.009 radians. If the seeker were 5 NMI from the target, this would translate to a cross 
range displacement error in the derived target position of roughly 83 meters assuming the 
seeker measured the azimuth and elevation precisely. 

As mentioned in Chapter II, DGPS, theoretically, should provide positions to 
within a one meter accuracy. Figures 4.1 through 4.3 show the GPS position errors 
relative to a NED orientation. From these Figures one can conclude the maximum 
amplitude for the displacement error occurs at roughly 12 seconds into the simulation over 
the sample period. At this point the value for this amplitude is about 2.5 meters. These 
errors are roughly consistent with actual DGPS performance specifications. 
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Figure 4.1: X-component of the GPS displacement error 




Figure 4.2: Y-component of the GPS displacement error 
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Figure 4.3: Z-component of the GPS displacement error 



The errors in the INS values for pitch roll and yaw are depicted in Figures 4.4 
through 4.6. From Figures 4.4 through 4.6, the maximum error in the INS’s computation 
of the Euler angles is approximately +/- 0.4 degrees. This equates to a cross range 
displacement error in the derived target position of approximately 65 meters at a 5 NMI 
range. The errors in the Euler angle computation for this model are comparable to the 
errors found from an actual INS system. Specifically, NRL’s P-3 uses the LTN-72 which 
is an INS with an advertised maximum amplitude for the Euler angle measurement error of 
less than 0.5 degrees. [Ref. 8, p. 35] 
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Figure 4.4: INS pitch error of the P-3 




Figure 4.5: INS roll error of the P-3 
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Figure 4.6: INS yaw error of the P-3 



B. SEEKER POSITION ERRORS 

Figures 4.7 through 4.12 show the position errors of seeker #1 and seeker #2. 
Notice, that these curves follow the errors associated with GPS very closely with slight 
variations due to the INS errors. As mentioned earlier, these errors translate to the same 
order magnitude errors in the derived target positions. There are two contributions to the 
overall error in deriving the position of each seeker. The first is the contribution from 
DGPS. The uncertainty in the DGPS position at the GPS antenna is lever arm corrected to 
the seeker position. In other words, the position errors of DGPS translate directly into 
position errors of the seekers. 
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Additionally, any angular errors from INS contribute to the overall position error 
of each seeker. A 0.5 degree error in any one of the Euler angle computations results in a 
5 centimeter error in the position of each seeker for a lever arm of 5 meters. This means 
that the errors from INS are insignificant in deriving each seeker’s position. This does not 
mean that the Euler angles are themselves are insignificant. Without the lever arm 
correction, the errors in deriving the seeker’s position in a fixed inertial coordinate system 
would be substantially greater. 




Figure 4.7: Seeker #1 position error in NorthySouth direction 
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Figure 4.8: Seeker #1 position error in East/West direction 




Figure 4.9: Seeker #1 position error in Up/Down direction 
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gure 4.10: Seeker #2 position error in North/South direction 




Figure 4.11: Seeker #2 position error in East/West direction 
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Figure 4.12: Seeker #2 position error in Up/Down direction 



Figures 4.13 and 4.14 provide a comparison of the true and derived positions of 
seeker #1. In Figures 4.13 and 4.14 the initial displacement in the +y direction for the 
position of seeker #1 is due to the location of seeker #1 in body fixed coordinates relative 
to the center of the P-3. Additionally, the P-3 does not fly a straight line. The course 
variations are due to the “P-3 FLIGHT PROFILER” model. Because this model 
continuously tries to reduce the error between the ordered course and the actual course, the 
P-3 can not maintain a steady course but continuously hunts for the ordered course. 

Figures 4.15 and 4.16 provide a comparison of the true and derived positions of 
seeker #2. The initial displacement in the -y direction is due to the placement of seeker #2 
with in the body fixed coordinates. 
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Figure 4.13: True position of seeker #1 




Figure 4.14: INS/GPS derived position of seeker #1 
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Figure 4.15: True position of seeker #2 




Figure 4.16: INS/GPS derived position of seeker #2 
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C. SEEKER CONTRIBUTED ERRORS 



Figure 4.17 shows the error from the range resolution for the conical scan seeker 
model used in the simulation. Because the range resolution is modeled as 100 meters, an 
envelope of 50 meters is evident throughout the simulation. Notice that the range error 
does not decrease over the sample interval. This is due to the fact that the range resolution 
is a function of the seeker hardware and not the signal to noise ratio. The range resolution 
of this seeker model translates to down range errors in the derived target position as is seen 
in the following section. 



1 1 r 




_j i i i i 

10 20 30 40 50 60 



Time (second) 

Figure 4.17: Range error from conical scan seeker (COSRO) 

Figure 4.18 shows the precision in the measurement of azimuth and elevation 
angles from the conical scan radar. This precision is approximated as +/- 0.15 degrees 
over this data collection period. This precision in the measurement of the azimuth and 
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elevation remains relatively constant during this period because the aircraft is closing the 
target from 5 NMI. Because the signal to noise ratio diminishes the effects of the error 
slope and crossover loss in the conical scan seeker, at this range, the errors in the azimuth 
and elevation are dominated by the scintillation error. 

An angular precision of +/- 0.15 degrees is expected to result in cross-range errors 
on the order of 25 meters when the seeker to target separation is 5 NMI. 
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Figure 4.18: Azimuth and elevation error from conical scan seeker (COSRO) 

Figure 4.19 shows the error in the range measurement for the monopulse seeker or 
seeker #2. Here the range resolution is 20 meters. As expected, the magnitude of the 
range error is +/- 10 meters. Notice the crossover rate for the range error is 5 times faster 
than for the crossover rate of the conical scan range error. This results from having a 
range resolution that is 5 times smaller in the case of the monopulse seeker. Again, the 
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range resolution of the monopulse seeker results in down range errors in the derived target 
position of roughly the same order magnitude. 




10 20 30 40 50 60 

Time (second) 

Figure 4.19: Range error from monopulse seeker 



Figure 4.20 shows the azimuth and elevation measurement errors for the 
monopulse seeker. Because the monopulse seeker is not affected by scintillation errors, an 
increase in signal to noise ratio directly results in an increase in the precision of the 
measurement. As the range between the seeker and target decreases, the azimuth and 
elevation measurement errors decrease from +/- 0.05 degrees to a negligible amount. At 
0.05 degrees the cross range error in the derived target position would be roughly 8 meters 
with a seeker to target separation of 5 NMI. 
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Figure 4.20: Azimuth and elevation error from monopulse seeker 



D. DERIVED TARGET POSITION ERRORS 

This section shows the true and derived target positions for each seeker type in the 
tangent plane coordinate system for several different types of target motion. In addition, 
this section shows the position errors that result from the DSTSPI algorithms for each 
case. 

Figures 4.21 through 4.29 were generated with the target dead in the water. 
Figures 4.30 through 4.38 were generated with the target on a course of 270° true. Figures 
4.39 through 4.47 were generated with the target on a course of 090° true. Figures 4.48 
through 4.56 were generated with the target on a course of true North. Figures 4.57 
through 4.65 were generated with the target on a course of 180° true. Figures 4.66 through 
4.74 were generated with the target making a slow turn to port from 000° to 270° true. 
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Figures 4.75 through 4.83 were generated with the target making a slow turn to starboard 
from 000° to 090° true. Figures 4.84 through 4.93 were generated with the target making 
a slow turn to starboard from 180° to 270° true. Figures 4.94 through 4.101 were 
generated with the target making a slow turn to port from 180° to 090° true. 

Table 4.1 is provided as a summary of the errors from the DSTSPI algorithms for 



seeker #1 and seeker #2 for each of the nine types of target motion. 





SEEKER #1 (COSRO) 


SEEKER #2 (MONOPULSE 


TARGET 


X 


Y 


Z 


X 


Y 


Z 


MOTION 


maximum 


maximum 


maximum 


maximum 


maximum 


maximum 


TYPE 


error 


error 


error 


error 


error 


error 




(meters) 


(meters) 


(meters) 


(meters) 


(meters) 


(meters) 


Dead-In-Water 


25 


42 


38 


7.5 


45 


40 


270° steady 


25 


48 


38 


7.5 


46 


40 


090° steady 


25 


40 


38 


8.1 


38 


40 


000° steady 


22 


43 


39 


11.5 


45 


41 


1 80° steady 


36 


43 


36 


9 


45 


26 


000° to 270° 


25 


42 


39 


10.5 


45 


40 


000° to 090° 


25 


42 


39 


10.5 


45 


40 


180° to 270° 


32 


42 


41 


10.5 


45 


40 


180° to 090° 


32 


42 


41 


10.5 


45 


40 


mean value 


27.44 


42.66 


38.77 


9.51 


44.33 


38.56 


std deviation 


4.66 


2.18 


1.56 


1.51 


2.39 


4.72 


variance 


21.77 


4.75 


2.44 


2.28 


5.75 


22.28 


maximum 


36 


48 


41 


11.5 


46 


41 



Table 4.1: Summary of DSTSPI derived target position errors 

From Table 4.1 several conclusions are made. The errors between the derived 
target positions from both seekers and the true target position are essentially constant over 
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a wide range of target motions. The errors are greatest in the cross range direction. The 
errors from the COSRO seeker are larger than for the monopulse seeker. This is especially 
true when the seeker-to-target separation is small. 

Assume that the errors that from the DSTSPI algorithms are roughly the same as 
the maximum values seen in Table 4.1. If this assumption is used, then the area of 
uncertainty in deriving a target’s position is roughly a sphere of radius 40 meters around 
the true target position. If the target’s speed is 22 kts or 11.3178 meters/sec, then the 
elapsed time before any decoys are detected after launch is about 3.5 seconds. This 
latency in detecting the decoy is defined as the time resolution, 8t, of the DSTSPI 
algorithms. The time resolution, 5t, decreases as the area of uncertainty of the target’s 
position decreases. From the figures presented in this section, the general trend of the area 
of uncertainty of the target’s position is that it decreases with a decreasing seeker-to-target 
separation. 
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Figure 4.21: True target position - DIW 
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Figure 4.22: Simulator #1 DSTSPI derived target position - DIW 
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Figure 4.23: Simulator #2 DSTSPI derived target position - DIW 
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Figure 4.24: Simulator #1 derived target position error in the North/South direction - 

DIW 








1 

TA 


i 

RGET 


i 

DEAD 


1 1 

N THI 


1 

i wat: 


ER 


1 1 


















































✓ 

s 


TRUE 

POSIT 


TARGET 

ION 










S 


s 


? 

KTS 




i - 








EED = 0 


























































1 


i i 




i 


i 


i 


i 


0.0 


NMI = 18.52 m 


j 


i i 



86 




Figure 4.25: Simulator #1 derived target position error in the East/West direction 
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Figure 4.27: Simulator #2 derived target position error in the North/South direction 

- DIW 




Figure 4.28: Simulator #2 derived target position error in the East/West direction - 

DIW 
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Simulator #2 derived target position error in the Up/Down direction - 

DIW 
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Figure 4.30: True target position - CSE 270° 
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Figure 4.31 : Simulator #1 DSTSPI derived target position - CSE 270° 




Figure 4.32: Simulator #2 DSTSPI derived target position - CSE 270° 
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Figure 4.33: Simulator #1 derived target position error in the North/South direction 

CSE 270° 




Figure 4.34: Simulator #1 derived target position error in the East/West direction 

CSE 270° 
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Figure 4.35: Simulator #1 derived target position error in the Up/Down direction - 

CSE 270° 
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Figure 4.36: Simulator #2 derived target position error in the North/South direction 

- CSE 270° 
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Figure 4.37: Simulator #2 derived target position error in the East/West direction 

CSE 270° 




Figure 4.38: Simulator #2 derived target position error in the Up/Down direction 

CSE 270° 
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Figure 4.39: True target position - CSE 090° 
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Figure 4.40: Simulator #1 DSTSPI derived target position - CSE 090° 
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Figure 4.41: Simulator #2 DSTSPI derived target position - CSE 090° 




Figure 4.42: Simulator #1 derived target position error in the North/South direction 

CSE 090° 
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Figure 4.43: Simulator #1 derived target position error in the East/West direction 

CSE 090° 




Figure 4.44: Simulator #1 derived target position error in the Up/Down direction 

CSE 090° 
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Figure 4.45: Simulator #2 derived target position error in the North/South direction 

- CSE 090° 





Figure 4.46: Simulator #2 derived target position error in the East/West direction 

CSE 090° 
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Figure 4.47: Simulator #2 derived target position error in the Up/Down direction 

CSE 090° 
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Figure 4,48: True target position - CSE 000° 
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Figure 4.49: Simulator #1 DSTSPI derived target position - CSE 000° 




Figure 4.50: Simulator #2 DSTSPI derived target position - CSE 000° 



99 




Figure 4.51: Simulator #lderived target position error in the North/South direction 

CSE 000° 




Figure 4.52: Simulator #1 derived target position error in the East/West direction - 

CSE 000° 
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Figure 4.53: Simulator #1 derived target position error in the Up/Down direction - 

CSE 000° 
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Figure 4.54: Simulator #2 derived target position error in the North/South direction 

- CSE 000° 
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Figure 4.55: Simulator #2 derived target position error in the East/West direction 

CSE 000° 




Figure 4.56: Simulator #2 derived target position error in the Up/Down direction 

CSE 000° 
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X AXIS = NORTH (NMI) x AXI S = NORTH (NMI) 
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Figure 4.57: True target position - CSE 180° 
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Figure 4.58: Simulator #1 DSTSPI derived target position - CSE 180° 
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Figure 4.59: Simulator #2 DSTSPI derived target position - CSE 180° 
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Figure 4.60: Simulator #1 derived target position error in the North/South direction 

CSE 180° 
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Figure 4.61 : Simulator #1 derived target position error in the East/West direction 

CSE 180° 




Figure 4.62: Simulator #1 derived target position error in the Up/Down direction 

CSE 180° 
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Figure 4.63: Simulator #2 derived target position error in the North/South direction 

- CSE 180° 




Figure 4.64: Simulator #2 derived target position error in the East/West direction - 

CSE 180° 
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X AXIS = NORTH (NMI) Z COMPONENT OF ERROR (meters) 




Simulator #2 derived target position error in the Up/Down direction - 

CSE 180° 




Figure 4.66: True target position - CSE 000° To 270° 
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Figure 4.67: Simulator #1 DSTSPI derived target position - CSE 000° To 270° 
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Figure 4.68: Simulator #2 DSTSPI derived target position - CSE 000° To 270° 
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Figure 4.69: Simulator #1 derived target position error in the North/South direction 

CSE 000° To 270° 




Figure 4.70: Simulator #1 derived target position error in the East/West direction 

CSE 000° To 270° 
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Figure 4.71: Simulator #1 derived target position error in the Up/Down direction 

CSE 000° To 270° 
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Figure 4.72: Simulator #2 derived target position error in the North/South direction 

- CSE 000° To 270° 
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Figure 4.73: Simulator #2 derived target position error in the East/West direction 

CSE 000° To 270° 




Figure 4.74: Simulator #2 derived target position error in the Up/Down direction 

CSE 000° To 270° 
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Figure 4.75: True target position - CSE 000° To 090° 
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Figure 4.77: Simulator #2 DSTSPI derived target position CSE 000° To - 090° 




Figure 4.78: Simulator #lderived target position error in the North/South direction - 

CSE 000° To 090° 
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Figure 4.79: Simulator #1 derived target position error in the East/West direction 

CSE 000° To 090° 




Figure 4.80: Simulator #1 derived target position error in the Up/Down direction 

CSE 000° To 090° 
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Figure 4.81: Simulator #2 derived target position error in the North/South direction 

- CSE 000° To 090° 




Figure 4.82: Simulator #2 derived target position error in the East/West direction 

CSE 000° To 090° 
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Figure 4.83: Simulator #2 derived target position error in the Up/Down direction - 

CSE 000° To 090° 




116 




Figure 4.85: Simulator #1 DSTSPI derived target position - CSE 180° To 270° 
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Figure 4.87: Simulator #1 derived target position error in the North/South direction 

CSE 180° To 270° 




Figure 4.88: Simulator #1 derived target position error in the East/West direction - 

CSE 180° To 270° 
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Figure 4.89: Simulator #1 derived target position error in the Up/Down direction - 

CSE 180° To 270° 




Figure 4.90: Simulator #2 derived target position error in the North/South direction 

-CSE 180° To 270° 
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Figure 4.91 : Simulator #2 derived target position error in the East/West direction 

CSE 180° To 270° 




Figure 4.92: Simulator #2 derived target position error in the Up/Down direction 

CSE 180° To 270° 
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X AXIS = NORTH (NMI) X AXIS = NORTH (NM1) 




Figure 4.93: True target position - CSE 180° To 090° 




Figure 4.94: Simulator #1 DSTSPI derived target position - CSE 180° To 090° 
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Figure 4.95: Simulator #2 DSTSPI derived target position - CSE 180° To 090° 




Figure 4.96: Simulator #1 derived target position error in the North/South direction - 

CSE 180° To 090° 
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Figure 4.97: Simulator #1 derived target position error in the East/West direction 

CSE 180° To 090° 




Figure 4.98: Simulator #1 derived target position error in the Up/Down direction 

CSE 180° To 090° 
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Figure 4.99: Simulator #2 derived target position error in the North/South direction 

- CSE 180° To 090° 




Figure 4.100: Simulator #2 derived target position error in the East/West direction - 

CSE 180° To 090° 
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Figure 4.101 : Simulator #2 derived target position error in the Up/Down direction - 

CSE 180° To 090° 
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V. HARDWARE IMPLEMENTATION 



As discussed in Chapter III, the DSTSPI algorithms require specific inputs. These 



inputs are a GPS position relative to the location of the GPS antenna, the pitch, roll and 



yaw of the aircraft that is given by the INS, and the azimuth, elevation and range between 



the seeker and the target. These inputs are shown in Figure 5.1. Of the three inputs, only 



the GPS position presents any real problems. 
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Figure 5.1: Hardware implementation of the DSTSPI algorithms 

The algorithms presented in Chapter III assumed all fixed inertial positions are 



given in tangent plane coordinates. This is done as a matter of convenience. To continue 



to use the tangent plane coordinates in the hardware implementation, a fixed position 



located in geodetic coordinates must be defined somewhere in the ASW/EW engagement 



test scenario. This fixed point of reference could be located at the center of the operating 
area that is allocated prior to the test. This point is taken from the existing navigational 
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chart, the coordinates of which are given in the geodetic sense of Latitude and Longitude. 
The actual location is arbitrary. 



To locate the P-3 in tangent plane coordinates we must take the position of the 
origin of the tangent plane coordinate system given in geodetic coordinates and perform a 
coordinate transformation from geodetic to ECEF coordinates. This position now 
becomes the origin of the tangent plane coordinate system given in ECEF coordinates. 

The geodetic to ECEF coordinate transformation is easily computed. First, three 
quantities must be defined. The first of these quantities is the flattening factor, f, which 
represents the relative flatness of the reference ellipsoid. The mathematical definition for f 
is defined as [Ref. 1, p. 118-119] 



where a and b are the semi -major and semi-minor axes of the reference ellipsoid, 
respectively. The next quantity that must be defined is the eccentricity, 8. The 
eccentricity of the reference ellipsoid is defined as 



The last quantity that must be defined is the normal, N, which represents the length of the 
ellipsoidal normal from the ellipsoidal surface to the intersection of the polar axis (Z ECEF ). 
Mathematically, N is defined as 



where <j) is the geodetic latitude. Now, by using these three quantities, the geodetic to 
ECEF transformation is defined by the following equations 



(5-2) 



N- 



a 



(5-3) 
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X ecef = (N + h) cos <|> cos X 


(5-4) 


Y ecef = (N + h) cos <|> sin X 


(5-5) 


z ecef =[N(l-E 2 ) + h]sin<t> 


(5-6) 



where h is the geodetic height or height above reference ellipsoid and X is the geodetic 
longitude. Equations (5-4) through (5-6) are incorporated into the MATLAB function 
“LL2ECEF.M” found in Appendix A. 

Now, the next step is to find the position of the P-3 relative to the tangent plane 
coordinates. To do this we take the GPS position given in ECEF coordinate and subtract it 
from the origin of the tangent plane given in ECEF coordinates. This results in a position 
vector that has the correct orientation, NED. 

Another option is to use a DGPS reference station. Here a portable DGPS 
reference station is located as near as possible to the testing location (this location must be 
on land). This station provides the necessary differential corrections to the local receiver. 
The advantage of using a portable DGPS reference station is that they usually provide 
better corrections than the Coast Guard differential stations. Another advantage is that 
these remote stations have the ability to record the GPS data. This would allow for further 
reductions in the errors from GPS. [Ref. 3] A typical system is found in Appendix B. The 
problem is that these reference stations must be setup and in operation 24 hours prior to 
the test scenario. Another problem is the cost. A typical reference station costs 
approximately twenty thousand dollars. 

Since the algorithms presented in Chapter III are not intended for real time use, an 
actual hardware implementation only involves recording the three inputs for post-test 
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analysis. The difficulty in doing this is to ensure the data is saved in a format that is 
compatible with a particular post analysis software. To make the process easier, several 
GPS receiver manufacturers provide software packages that provide a post analysis 
capability. To this end it behooves the user to pick a GPS system that comes with a 
complete software package. 

Another potential hazard in the recording of the data for post-test analysis is to 
ensure that all the recorded times are properly synchronized. Any clock differences 
translate directly into position errors. This is especially true for the time-tagged data of 
the INS and the Simulators since the angular rates can be quite large. 

Perhaps the most difficult problem is the determination of the actual target position 
as a function of time given noisy data from the DSTSPI algorithms. The goal of the post- 
test analysis is to provide sufficient accuracy in the determination of the target’s position 
so as to allow for the recognition of any launched decoys. To help attain this goal it is 
possible to take the raw data from the DSTSPI algorithms and in some way make an 
educated guess about the true target position. A first order approximation in obtaining the 
target’s true position might be a simple least-squares fit of the data, but because of the 
cross range errors which are a consequence of the INS and seeker errors, a simple least- 
squares fit on the data does not provide any useful information. This is especially true for 
a maneuvering target. A better method in determining the target’s true position given the 
noisy data might involve an iterative process by which the errors in all directions are 
minimized at the same rate. 
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VI. CONCLUDING REMARKS 



The effectiveness of the at sea Anti-Ship Missile/Electronic Warfare engagement 
test scenario is heavily dependent upon the results of this thesis. By knowing the target's 
position in a fixed inertial coordinate system, the post analysis problem of distinguishing 
between target and decoy is made easier. Any launched decoys during the test scenario 
show up as an instantaneous deviation in the target track. This allows the test engineers 
determine if the seeker broke lock on the target and tracked an offboard decoy (e.g., chaff). 
By improving the method by which the at sea ASM/EW test scenario evaluation is 
performed, the war-fighting capabilities of the U. S. Navy improve. 

Also, by knowing a target's position relative to a fixed inertial coordinate system, a 
target's position can be sent to any platform and be understood immediately without being 
specific to the platform’s Fire Control System. Target positions that are fixed in a fixed 
inertial coordinate system allow for ease of processing because the information is sensor 
independent. 

Considering the results obtained in this simulation, it is feasible to derive a target’s 
position in a fixed inertial coordinate system from the integration of INS, GPS and radar. 
The precision to which the target’s position is derived depends heavily upon the quality of 
the navigation equipment used in the integration process. 

With the advent of more and more sophisticated weaponry, it becomes increasingly 
more important to provide the best available means by which an evaluation of these 
weapons is made. Of particular significance is the foreign missile threat. It is here that a 
great deal of time and money is spent in determining the effectiveness of our shipboard 
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self defense systems against these threats. The validation of a shipboard self defense 
system against a missile threat usually involves modeling and simulation. 
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APPENDIX A 



MATLAB FILES: 



AZIMUTH.M 



function out = azimuthfthetah 
%ANDREW ROWE 

%THIS FUNCTION COMPUTES THE AZIMUTH ANGLE 

%RELATIVE TO TRUE NORTH FROM THE THETA GIVEN 

%FROM CARTESIAN TO SPHERICAL COORDINATE 

%TRANSFORMATION 

%1 0/1 0/96 

x = theta; 

y=450-x; 

if y >= 360 

y=y-360; 

end; 

out = y; 



B2NANT.M 



function out = b2nant(vecl: 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 
% ROLL=PHI PITCH=THETA YAW=PSI 

% pi IS THE POSITION VECTOR OF THE GPS ANTENNA IN BODY FIXED FRAME 

% COORDINATES RELATIVE TO THE CENTER OF THE AIRCRAFT 

% 10/10/96 

phi = vec(l); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos( theta) 0 sin(theta); 0 10; -sin(theta) 0 cos(theta)]; 

Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1]; 
pl=[-5;0;0]; 

out = inv(Tpsi*Ttheta*Tphi)*pl; 



B2NCEN.M 



function out = b2ncen(vecl: 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pi IS THE POSITION VECTOR OF THE CENTER OF THE AIRCRAFT RELATIVE TO THE 

% POSITION OF THE GPS ANTENNA 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(l); 

theta = vec(2); 
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psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos(theta) 0 sin(theta); 0 10; -sin( theta) 0 cos( theta)]; 
Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1]; 
pl=[5;0;0]; 

out = mv(Tpsi*Ttheta*Tphi)*pl; 



B2NP1.M 



function out = b2np lived: 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pi IS THE POSITION VECTOR OF THE AFT PART OF THE SIMULATOR #1 RELATIVE TO 

% THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(l); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos(theta) 0 sin(theta); 0 10; -sin( theta) 0 cos(theta)]; 

Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1]; 
pl=[-2;5;0]; 

out = inv(Tpsi*Ttheta*Tphi)*pl; 



B2NP12.M 



function out = b2nnl2('vec): 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pi IS THE POSITION VECTOR OF THE AFT PART OF THE SIMULATOR #2 RELATIVE TO 

% THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(l); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos( theta) 0 sin(theta); 0 10; -sin(theta) 0 cos(theta)]; 

Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1]; 
pl=[-2;-5;0]; 

out = inv(Tpsi*Ttheta*Tphi)*pl; 



B2NP2.M 



function out - b2np2(vecl: 

% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pi IS THE POSITION VECTOR OF THE FORWARD PART OF THE SIMULATOR #1 
% RELATIVE TO THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 
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% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(l); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos( theta) 0 sin(theta); 0 10; -sin(theta) 0 cos(theta)]; 
Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1]; 
pl=[2;5;0]; 

out = inv(Tpsi*Ttheta*Tphi)*pl; 



B2NP22.M 



function out = b2np22(vecL 
% ANDREW ROWE 

% THIS FUNCTION CONVERTS A VECTOR GIVEN IN THE BODY FIXED FRAME COORDINATE 
% SYSTEM INTO THE NAVIGATION FRAME 

% pi IS THE POSITION VECTOR OF THE FORWARD PART OF THE SIMULATOR #2 

% RELATIVE TO THE CENTER OF THE AIRCRAFT IN BODY COORDINATES 

% ROLL=PHI PITCH=THETA YAW=PSI 

% 10/10/96 

phi = vec(l); 

theta = vec(2); 

psi = vec(3); 

Tphi = [1 0 0; 0 cos(phi) -sin(phi); 0 sin(phi) cos(phi)]; 

Ttheta = [cos(theta) 0 sin(theta); 0 10; -sin(theta) 0 cos(theta)]; 

Tpsi = [ cos(psi) -sin(psi) 0; sin(psi) cos(psi) 0;0 0 1]; 
pl=[2;-5;0]; 

out = inv(Tpsi*Ttheta*Tphi)*pl; 



CONTHETA.M 



function out = conthetaf component); 

% ANDREW ROWE 

% THIS FUNCTION COMPUTES THE ANGULAR PRECISION FOR A CONICAL SCAN RADAR 
% SEEKER. THIS IS USED IN THE SEEKER MODEL 
% 10/10/96 



R = component 1); %range in meters 

Pt = 250e3; %power in watts 250 Kw 

sigma = 3000; %radar cross section square meters 

f = 9e9; %radar frequency of 9 Ghz 

alpha = .055; %attenuation dB/km page279 Barton 

Gain = 1995.26; % antenna gain 33 dB 

c = 2.9979e8; %speed of light 

lambda = c/f; %wavelenght in meters 

B = 10e6; %noise bandwidth 

F = 10; %noise figure 

kT = 4. 14e-21 ; %joules at 300 K 

bw = .0785398; %3dB transmitter beam width (uniform illumination) 4.5 degrees 

SN = Pt*Gain*Gain*sigma* 1 0 A (-0.2*alpha*R/l 000)* lambda* lambda/4/pi/kT/B/F/4/4/pi/pi/R/R/R/R; 

sigma_s=0.0125*bw; %scintillation error at optimized squint angle 

%Barton page 390 

SN = 100*SN; % 100 pulse integration improvement 
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thetasig = bw/(SN) A 0.5 + sigma_s; 



%radians 



THETAD.M 



function thetasig = thetad(component); 

% ANDREW ROWE 

% THIS FUNCTION COMPUTES THE ANGULAR PRECISION FOR A MONOPULSE 
% SEEKER. THIS IS USED IN THE SEEKER MODEL 
% 10/10/96 



R = component 1); 
Pt = 30e3; 
sigma = 3000; 
f = 17e9; 
alpha = .055; 

Gain = 218.77; 
c = 2.9979e8; 
lambda = c/f; 

B = 10e6; 

F = 8.9; 
kT = 4.14e-21; 
bw = . 13962; 



%range in meters 

%power in watts 30 Kw 

%radar cross section square meters 

%radar frequency of 17 Ghz 

%attenuation dB/km page279 Barton 

% antenna gain 23.4 dB 

%speed of light 

%wavelenght in meters 

%noise bandwidth 

%noise figure 



%joules at 300 K 

%3dB transmitter beam width (uniform illumination) 8 degrees 
SN = Pt*Gain*Gain* sigma* 10 A (-0.2*alpha*R/1000)*lambda*lambda/4/pi/kT/B/F/4/4/pi/pi/R/R/R/R; 
km = 1 .9; %error slope for monopulse radar dual horn 

thetasig = bw/km/(100*2*SN) A 0.5; %radians 



HEAD.M 



function out = head(direc): 

% ANDREW ROWE 

% THIS FUNCTION DETERMINES THE HEADING BASED ON NOT EXCEEDING 180 DEGREE 
% TURN DURING AIRCRAFT MANEUVERS 
% 10/10/96 
x = direc(l); 
y = direc(2); 

if (x > y) & ((x-y) < pi), 
direction = 1; 

elseif (x < y) & ((y-x) < pi), 
direction = -l; 

elseif (x > y) & ((x-y) >= pi), 
direction = - 1 ; 

elseif (x < y) & ((y-x) >= pi), 
direction = 1; 

elseif x=y, 

direction = 1; 

end; 

out = direction; 



OVER2PLM 



function out = over2pi(direc); 

% ANDREW ROWE 

% THIS FUNCTION ENSURES THAT THE HEADING FALLS BETWEEN 0 AND 2PI 
% 10/10/96 
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LL2ECEF.M 



x = direc; 
if x >= 2*pi 

x=x-2*pi; 

end; 
if x < 0, 

x=x+2*pi; 

end; 
out = x; 



function ecef=l 12ecef(x) 

% converts from geodetic latitude, longitude, elevation to % earth-centered earth-fixed Cartesian 
coordinates 

phi = x(l); lambda = x(2);h = x(3); 

% define semi-major and semi-minor earth axes a=6378137;b=6356000; 

% define auxiliary quantities f,e, and N 
f=(a-b)/a; e = F(2-f); N = a/sqrt(l-(e*sin(phi))-2); 

% convert to Cartesian 

ecef = [(N + h)*cos(phi) * cos(lambda) 

(N + h)*cos(phi) * sin(lambda) 

(N*(l-e A 2)+h) * sin(phi)]; 
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APPENDIX B 



The following Pages include manufacturer data-sheets from the DEL NORTE 
Technology, Inc. 

DEL NORTE Technology. Inc . e-mail: dnti@delnorte.com 
1) FLYING FLAGMAN® Global Positioning System Measuring Unit 

FEATURES/BENEFITS : 

• 12-Channel GPS Receiver 

• Flow Rate/Temp/Humidity Monitoring 

• 10 Home Points 

• 100 Field Memory 

• Polygon Field Entry 

• Racetrack, Parallel Line and Squeeze to Middle Functions 

• 10 Position/Sec Updates 

• Ground Speed in MPH, m/sec or knots 

• Serial Port for control/logging 

• Break Point Memory 

• Resume to Point 

• Automatic Reverse Direction Sensing 

• Skipped Line Detection 

• Selection of feet/meters 

• Built-in Clock 

• Data/Time Tag 

• Remote Functions on Flight Stick 

• On-Ground Self-test 
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The new DEL NORTE Technology GPS Measuring Unit (GMU) provides a 
rugged, waterproof, and versatile GPS Receiver with multiple function capabilities and a 
range of options Set up is accomplished by using a menu driven Control and Display Unit 
(CDU). The Del Norte GMU can serve as either a Reference or Mobile GPS Unit. As a 
Reference Unit, the GMU, with its Establish Position capability, is able to derive site 
location to within one meter in a stand-alone mode over twenty-four hours and even better 
if left for forty-eight hours. 

After setup the CDU may be disconnected to leave the single rugged, waterproof 
unit on site. By adding a 'second processor' option to the GMU, the Reference Station can 
also act as a differential monitor. As a mobile, the GMU may be interfaced to any of the 
user's data processing systems. For stand-alone operations the addition of the VGA board 
and optional second processor allow the GMU to drive a color display and incorporate 
third party software. 

SPECIFICATIONS : 

PERFORMANCE : 

• Differential Dynamic: 1-3 meters, 1 meter CEP 

• Update Rate : 0.1 seconds 

SYSTEM ENVIRONMENT : 

• Temperature: 0-55 degrees Celsius (32 - 130 degrees Fahrenheit) 

• Vibration: Tested to MIL-STD-810-C 
PHYSICAL/POWER CHARACTERISTICS : 

• System Weight: 24 lbs (10.88 Kg) • Super Trac Light Bar 

• CDU Dimensions • W 15.0 in (38.1 cm) 

• W 4.2 in (10.6 cm) • H 2.0 in (5.1 cm) 

• H 7.7 in (19.5 cm) • D 5.5 in (14.0 cm) 
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• D 1.3 in (3.3 cm) 

• GMU Dimensions 

• W 8 in (20.3 cm) 

• H 6.8 in (17.3 cm) 

• D 15 in (38.1 cm) 

SYSTEM SUPPLIED WITH 

• GMU Mainframe 

• Control/Display Unit 

• Super Trac Light Bar 

• Remote Step Switch 

• GPS Antenna 
OPTIONS : 

• Differential Data Sources 

• Local Reference Tower 

• Portable Reference Station 

• FM Side-Band 

• Government MF Transmissions 

• Satellite Data Links 

• Data Logger and Post-Processing/Plotter Software 

• Map Trac Moving Map Display 
SUPER TRAC LIGHT BAR : 

• Pilot Programmable Displays • 

• Digital HSI 

• Line/Pass Number 

• Distance Off-Track 

• Distance-to-go • 

• Course Deviation 

• Area Sprayed 

• Custom Settings 

• Ground Speed 

2) MODEL 2012 PORTABLE REFERENCE STATION 



Left/Right Guidance 

• Visible in Direct Sunlight 

• Dimmer Control for Low Light/Night Flying 

• Programmable Sensitivity 
Course Deviation Indicators 

• Programmable to nearest 0.1 of a degree 



• Power Input : 11-32 VDC, 30 Watts 



• Differential Rx + Antenna 

• Cables 

• Shock Mounts 

• Installation and Operation Manuals 
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FEATURES/BENEFITS: 



• Dedicated GPS Reference Station • 12-Channel GPS Receiver and Radio/Modem 

• Lightweight & Portable which can be in one box 

• Low-Cost • 'Self-Locating' 

The Del Norte 2012 is the low-cost portable reference station solution. To 

compliment the 2012, a small, lightweight, and inexpensive UHF Mobile Data Link 
Transmitter and Receiver is available to feed Differential GPS Corrections to your mobile 
GPS receiver. Del Norte has trained engineers and technicians available for field 
deployment anywhere in the world to provide assistance with first-time installations or 
emergency needs 

SPECIFICATIONS : 

PHYSICAL & POWER CHARACTERISTICS : 

• Weight 5 Kg • 

• Dimensions 

W 30 cm • 

H 9 cm 
D 34 cm 

SUPPLIED WITH 

• 12-Channel LI Receiver • 

• UHF Transmitter and Modem • 

• Patch GPS Antenna w/Ground Plane • 

• Broad Band Omni Antenna 

• 

SYSTEM OPTIONS : 

• 21 5UL Mobile Data Link Receiver 
RADIO SPECIFICATIONS : 

• Frequency 

458.xxx MHz, 0.5 Watt (UK deregulated), options include other UHF, VHF, and 
HF frequencies on request 

• Transfer Times 



Operating Temperature 
0-50 Degrees Celsius 

Power 

11 - 14.5 vDC 

Necessary Cables and Connectors 
Operator's Manual 

Comprehensive Product and Firmware 
Maintenance for One Year 
Training Available 



144 



W/RTCM Format 

Typical w/4 satellites, 0.4 seconds 

Typical w/12 satellites, 1.0 seconds 

• Baud Rates 
300 to 9600 
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